Effective low-energy theory for correlated carbon nanotubes 
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The low-energy theory for single-wall carbon nanotubes including Coulomb interactions is derived and analyzed. It describes 
two fermion chains without interchain hopping but coupled in a specific way by the interaction. The strong-coupling properties 
are studied by bosonization, and consequences for experiments on single armchair nanotubes are discussed. 
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The recent discovery of carbon nanotubes (T) has 
sparked a tremendous amount of activity. Nanotubes 
are nanoscale particles obtained by wrapping a single 
layer of graphite into a cylinder. The electronic prop- 
erties of a (n, m) tube are determined by the integer in- 
dices < to < n of the wrapping superlattice vector, 
and depending on the choice of m and 77, the tube is 
either a metal, a narrow-gap semiconductor, or an in- 
sulator Experimentally, nanotubes can be produced 
using the carbon-arc technique Q or by laser ablation ||] 
of Co- or Ni-doped graphite targets. The latter method 
yields metallic (77, 77) "armchair" nanotubes with 77 = 10 
in rather large quantities, usually deposited in triangular- 
packed ropes. Remarkably, the first transport measure- 
ments on a single 3/im long (10, 10) nanotube have been 
reported recently j|. 

Carbon nanotubes are perfect experimental realiza- 
tions of one-dimensional (ID) conductors. Interacting 
ID electrons usually exhibit Luttinger liquid rather than 
Fermi liquid behavior characterized by, e.g., the absence 
of Landau quasi-particles, spin-charge separation, sup- 
pression of the electron tunneling density of states, and 
interaction-dependent power laws for transport quan- 
tities [|. So far non-Fermi liquid behavior has been 
masked by charging effects due to large contact resis- 
tances between the nanotube and the attached leads j| . 
Nevertheless, future experiments are expected to reveal 
the anomalous conductance laws and related phenomena 
discussed here at higher temperature, suitable gate volt- 
ages, or smaller contact resistances. 

The low-energy theory for uncorrelated nanotubes has 
been given by Kane and Mele |(|. Here we extend their 
approach and incorporate Coulomb interactions among 
the electrons, focusing on (77, 77) tubes where interaction 
effects are most pronounced. Previously, this problem 
has only been studied by the perturbative renormaliza- 
tion group (RG) using a weak short-range (Hubbard) in- 
teraction |7|. In this Letter, we discuss arbitrary inter- 
action potentials and, in particular, the strong-coupling 
regime which emerges at low temperatures. This is of 
importance as the experiments described in Ref. [Q are 



characterized by long-ranged (unscreened) Coulomb in- 
teractions. The low-energy theory found here is equiv- 
alent to two spin-i fermion chains coupled in a rather 
special way by the interactions but without interchain 
single-particle hopping. In that respect, our theory dif- 
fers from the standard two-chain problem [|]|| ■ 

The remarkable electronic properties of carbon nan- 
otubes are due to the special bandstructure of the it 
electrons in graphite The Fermi surface consists 

of two distinct Fermi points aK with K = (47r/3a, 0) 
and a = ±, see Fig. [y. Here the cc-axis points along 
the tube direction and the circumferential variable is 
< y < 2nR, where R = \/3na/2w is the tube ra- 
dius. The lattice constant is a = 2.46A. Since the ba- 
sis of graphite contains two carbon atoms, there are 
two sublattices p = ± shifted by the vector d = (0,d) 
[where d — a/y/3], and hence two degenerate Bloch states 
<Ppa(x,y) at each Fermi point aK. We follow Ref. Jl|] 
and choose these states separately on each sublattice such 
that they vanish on the other (note that K y = 0), 



(p pa (x,y) cp- pa >(x,y) = , ip pa {x,y) 



e~}q>(—iaK x x) 



At low energy scales, the electron operator for spin a = ± 
can be expanded in terms of the Bloch waves |(]Jll]] , 



* CT (X, y) = ^ Vpaix, y) ^paa(x) . 



(1) 



Quantization of transverse motion gives the wavefunc- 
tion x(y) = exp(iMy) with the "mass" M specified in 
Ref. |(|. Therefore \ = 1 for armchair tubes (M = 0), 
see Eq. ([!]). Excitation of other transversal bands costs 
the energy « 10 eV/77, and hence a ID situation arises. 
The conditions for the low-energy regime are met even 
at room temperature for (10, 10) nanotubes. Neglecting 
Coulomb interactions, the Hamiltonian is ]^JTl[| 



Ho = -v^p / dx ipl a(T dxip- po 



(2) 



where v ~ 6 x 10 6 m/sec is the Fermi velocity. Switching 
from the sublattice (p = ±) description to the right- and 
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left- movers (r = ±) indicated in Fig. [l] implies a massless 
ID Dirac Hamiltonian. 
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FIG. 1. Low-energy bandstructure of graphite. Fermi 
points are labeled by a = ±, and sublattices p = ± combine 
to build right- and left-movers (r = ±). 

Let us now examine Coulomb interactions mediated by 
a (possibly screened) potential U(x—x',y — y'). Electrons 
trapped in nonpropagating orbitals can be incorporated 
in terms of a dielectric constant k, but free charge carriers 
in nearby gates could lead to an effectively short-ranged 
potential. For the experiments of Ref. one has an 
unscreened Coulomb interaction, and the dielectric con- 
stant can be estimated astcw 1.4 @. From Eq. (§) the 
interaction contribution reads 
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* = \ £ E J dxdx ' v {U x - x ') ® 

pp' aa' aia.20-30.4 
X ^l ai a{x)^ a ^{x')-l}} pla3a ,{x)^ pai<J {x) 

with the ID interaction potentials 



V { Z } (x-x>) 



2nR 



dydy' v* pai (x, y)<p* p > a2 (x',y') (4) 



x U(x-x',y-y' +pd5 p -p>) <p p > , 3 (x',y')<p pai (x,y) . 

Here inter-sublattice interactions involve the shift vec- 
tor d = (0,d). It is natural to distinguish three pro- 
cesses associated with the Fermi points a = ±. First, 
we have "forward scattering" (oFS) where u\ = an and 
ot% = «3- Second, we have "backscattering" (aBS) with 
ai = —ai = 013 = — u±. Finally, at half-filling there 
is an additional "Umklapp" process (aUS) characterized 
by a\ — ct2 = — «3 = — «4. These processes are different 
from the conventional ones || since they do not necessar- 
ily mix right- and left-moving branches but rather involve 
different Fermi points, see Fig. |l|. 

Let us start with oFS. We first define 



V (x) 



2-kR 



dy 



2 * R dy' 



o 2ttR J 2ttR 



U(x,y-y'), (5) 



whence from Eq. (^) the {ai}-independent ID potential 
reads V^p S (x) — Vo(x) + 5 p .-p<8V p {x) with 



f 2vR dvdv' 

SV p = I T^h M*, y-y'+pd)-U{x,y- y')] 



(2itR) 2 



The correction SV P measures the difference between intra- 
and inter-sublattice interactions. Because of the period- 
icity of the y-integrals, expanding in powers of d implies 
8Vp{x) = 0. Since the potential Vq does not discriminate 
among sublattices, the resulting oFS interaction couples 
only the total ID charge densities, 



H 



(o) 

qFS 



With p{x) = J^paa^paa^paa 

Coulomb interaction (ao — a), 

U(x,y) = 



dxdx' p(x)Vo(x — x')p(x') , (6) 
For an unscreened 



n^a 2 + x 2 + AR 2 sin 2 (y/2i?) 
the potential (^) becomes 

2e 2 / 2R 



(7) 



V (x) 
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with the complete elliptic integral of the first kind K(z). 
The Fourier transform Vb(fc) for \kR\ <C 1 is reminiscent 
of a ID quantum wire B, 



V (k) = (e 2 /n) [2| ln(fci?)| + 7rln2] 



(8) 



For \x\ a, the above continuum argument leading to 
5V p (x) = safely applies. However, for |x| < a, an ad- 
ditional aFS term beyond Eq. (^) arises due to the hard 
core of the Coulomb interaction, 



H 



(i) 

aFS 



-f J dx ^paa^-pa'a'^-pa'a'^p 



(9) 



with f/a = -fe 2 /R. Evaluating SV p (x — 0) on the 
wrapped graphite lattice using Eq. (0) yields 



V3a 

27TKO0 



1 - 



y/i + a 2 /3a, 



0.1 



(10) 



In the language of a Hubbard model, we have f/a = 
U — V where U = e 2 /R is the on-site and V the nearest- 
neighbor Coulomb interaction. According to Eq. (|Tc|), 
this difference is small compared to U. 

Next we discuss aBS contributions. Since Eqs. @ and 
(J|) involve a rapidly oscillating factor exp[2iK x (x — x')], 
these are local processes which do not resolve sublattices, 



H a BS — 2 



dx 



£ ^paa^l'-CKT'^p'aa'^p-a* ■ (H) 

pp' OL<J<j' 
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Estimating the coupling constant b from Eq. (Q) for the 
unscreened interaction (Q) results in b « /, while for well- 
screened short-ranged interactions, b 3> /. Experimen- 
tally, by using additional gates, one can easily tune away 
from half-filling [Q . Therefore we disregard all Umklapp 
scattering effects in this paper. 

To study this effective low-energy model, we employ 
standard Abelian bosonization [pH. For that purpose 
we first switch to right- and left-movers (r = ±), 

= X! U w i>ra<j where U^OyU = a z , 

T" 

which then allow for straightforward bosonization, 




Phase fields for the total and relative (6 — ±) charge 
(j = c) and spin (J = s) channels obey the algebra 



j5 (x),e j , s >(x')} = -{i/2)5. 



of 



sgn(x 



Thus <pjs and 9js are dual fields. Density and current 
operators can be written in terms of the c+ field, 

p{x) = (2/V^)d x 9 c+ (x) , / = {2e/V^)d t 9 c+ (0) . 

The wavevector qp is related to deviations bp = Aqp/ir of 
the average electron density from half-filling and can be 
tuned by gates. Finally, the rj raiy are Majorana fermions 
ensuring anticommutation relations among different ran. 
Since the Hamiltonian contains only the combinations 
A±± = r/ raa -rj± r ± a<7 , these can be represented using 
standard Pauli matrices. Besides A ++ = 1, we have 
A_| = iaa x , A |_ = iracr z and A = —ir<j y . Us- 
ing Eq. (E2I), the bosonized expressions then read 



H = V - dx j: 



[(d x <t> jS ) 2 + Kr s 2 (d x 9 jS ) 2 ] (13) 



H aFS = ~ dxdx ' d x 9 c+ (x)V Q (x - x')d x ,9 c+ (x') (14) 



/ 



dx [— cos(v47r# c _) cos(v 47r # s _) 

- cos(\/47r 9 C -) cos(V4^0 s+ ) (15) 
+ cos(v / 47r6' s+ ) cos(\/47r# s _)] 

H a BS = 7 — 77 / dx[cos(V4^9 c -)cos(V47r9 s -) 
{■nay J 

+ cos(\/47r6» c _) cos(V4~7r </),,_) (16) 
+ cos(v / 47T0 s -) cos(v / 47T0 s -)] • 

Although bosonization of Eq. (|J) gives Kjg — 1 , interac- 
tions renormalize these parameters. In particular, if^pg 
can be incorporated into Hq by putting K c+ = K = 



l/\/l + 4Vo(fc — 0)/nv < 1, while for all other channels, 
Kjs = •v/l + / / ttv > 1. For the long-ranged interaction 
(0), one has K — const, j ^J\ ln[max(T, v/L)]\ || such that 
at zero temperature, K vanishes in a very long tube. Ve- 
locity renormalizations are ignored here as they do not 
change exponents. 

From the perturbative RG equations |L3|], it follows 
that around the Gaussian point / = b = 0, H^p S is 
irrelevant, while H a Bs scales to strong coupling, b — > 00. 
Let us first examine short-ranged interactions such that 
/ <C b. For / — 0, the relative charge and spin channels 
can be refermionized in terms of four Majorana fermions 
£jRtL)(x). Here j = 1,2 corresponds to (s— ), j = 3,4 
to (c— ), and R(L) denotes the right/left-moving part. 
Remarkably, refermionization shows that the last term 
in Eq. (fLq) vanishes identically, and we are left with 



h(c-,s-) = -^Yl / dx ,; s.//."'U> i;,/ <'u; ; /.; 



3=1 

-2b J dx (6a6l + UrUl)^il ■ (17) 

The Majorana fermion £2 stays massless, i.e., the (s— ) 
sector carries one massive and one massless branch. This 
behavior is due to the symmetric appearance of the dual 
fields 9 S - and </> s _ in Eq. (^|) which does not permit 
complete pinning due to the Heisenberg uncertainty re- 
lation The (c— ) sector is fully gapped, and we can 
put 9 C - = 0. The masses of the three massive Majorana 
fermions £1,3,4 at this strong-coupling point are approx- 
imately equal (sa nib) and can be estimated from mean- 
field theory as nib ~ lo c exp(— nv/y/2b), where u c = 7.4 
eV is the bandwidth of the 7r electrons 0. Correla- 
tion functions for the (s— ) sector follow from the cor- 
respondence between two Majorana fermions and the or- 
der/disorder operators of the 2D Ising model ||. In par- 
ticular, one has scaling dimension 77 = 1/8 for sin(y / 7r# s _) 
and sin(y / 7T0 s _) [their scaling dimension for b = is 
i] = 1/4], while correlation functions of cos(v / 7r6' s _) and 
cos(v / 7T0 s _) decay exponentially. In the (c— ) sector, cor- 
relations of cos(y / 7r0 c _) show long-range order while all 
other operators lead to exponential decay. 

While H^p S is irrelevant at / = b = 0, it becomes rel- 
evant at the new strong-coupling point which therefore 
describes only an intermediate fixed point. The term 
~ cos(V47r6' s +) cos(y/4ir9s-) in Eq. (|l^) stays marginal, 
but the two other terms become relevant with scaling 
dimension 77 = 1. Therefore £2 as well as the (s+) field 
acquire the mass rrif ~ (/ /b)mb- At the emerging strong- 
coupling fixed point, we have long-range order in the op- 
erators cos(^/^^9 s +), cos(y / 7r6' c _) and sin(y / 7r<?!>, s _), with 
exponential decay in all other operators except those of 
the (c+) sector. 

This analysis for a screened interaction potential pre- 
dicts that the exponents corresponding to the first (inter- 
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mediate) strong-coupling point should be observable on 
temperature scales m/ < T < nib with a crossover to a 
regime T < mt dominated by the true T = fixed point. 
For long-ranged interactions, we have nif « m&, and 
the intermediate fixed point and the associated crossover 
phenomenon are then absent [ff2|. 

With the strong-coupling solution discussed above, we 
can now examine temperature-dependent susceptibilities 
and other experimentally accessible quantities of inter- 
est. At temperatures T > mj, the dominant correla- 
tions come from the inter-sublattice charge-density wave 
(CDW) and spin-density wave (SDW) operators 



Ocdw 



• OsDW 



paa 



with (O(x)6(0)) ~ cos(2q F x)x- {K+ ^/ 2 for both ||. 
Due to a larger prefactor, SDW correlations will be more 
pronounced. For nif < T < nib, this decay is changed 
into the slower x ~(2K+3)/4 law ^ Both the CDW anc j SDW 

correlations decay exponentially at very low tempera- 
tures T < nif. However, there is also a iqp CDW com- 
ponent effectively coming from C p,w ll> which leads to 
a slow cos(4qp x)x~ 2K decay at T < nit,. This is in fact 
the dominant instability at temperatures T < nib and 
strong correlations, K < 1/2. Remarkably, it is not K x 
but <7f which determines the period of all the dominating 
correlations. Since aqp <C 1, one has pronounced ferro- 
magnetic correlations. This offers an explanation for the 
ferromagnetic tendencies observed in Ref. Q. 

Superconductivity (SC) has been predicted to be quite 
a robust feature of two-chain models ||7|-[j). In the nan- 
otube, the dominant SC correlations come from the intra- 
sublattice singlet SC pairing operator 



Osc 



such processes will not exert backscattering in the (c+) 
channel. Hence they do not affect the conductance. The 
latter two effects, however, mix r = ± modes and thus 
couple to the (c+) mode. Both lead to mass-term like 
perturbations B of the generic form 



dxM(x) 5] 



pa a 



Assuming a single static pointlike impurity center, in 
order M 2 the following temperature dependence of the 
corrections SG to Go is found. For T > nib, we obtain 
SG ~ T^'- 1 )/ 2 , which is turned into SG ~ T^-^' 4 at 
nif < T < nib- At even lower temperature, T < mf, 
from perturbation theory in M 2 , one would get 5G ~ 
T K ~ 2 . For extremely low temperatures, however, by in- 
voking the usual duality argument || , this is turned into 
G(T — > 0) ~ 7 1 - 2 + 4 / K j i. e . ; a vanishing conductance 
similar to the Luttinger liquid case. For long-ranged in- 
teractions, this leads to a pseudogap behavior at very low 
temperatures. 

To conclude, the effective low-energy theory for corre- 
lated single-wall nanotubes has been given. Our theory 
explains the ferromagnetic tendencies observed in recent 
experiments and predicts the temperature dependence of 
various susceptibilities and the conductance in the pres- 
ence of impurities. When additional gates are present, 
an interesting crossover related to screened interactions 
should be observable. 
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SRC of the United Kingdom. 



while triplet SC plays no role. For T > nib, correla- 
tions decay faster than x~ 2 , and for nif < T < nib, their 
a ,-(2/K+3)/4 decay is su bdominant to the CDW and SDW 

correlations. Finally, at very low temperatures, T < nif, 
we obtain (d S c(x)d S c{0)) ~ x-^ 2K . Therefore SC be- 
comes the dominant instability only at very low tem- 
peratures and for sufficiently short-ranged interactions 
(K > 1/2). For long-ranged interactions, superconduct- 
ing correlations are of no importance. 

The properties of the above strong-coupling points 
also determine conductance laws. In the absence of 
impurities, one has perfect conductance quantization, 
G = Go = 2e 2 /h, and experimentally measured resis- 
tances are contact resistances. There are nevertheless 
various sources for impurities, e.g., structural imperfec- 
tions, charge defects in the substrate, topological defects 
(dislocation pairs) |Ej, or twists M. Since the first two 
act like elastic potential scatterers coupling to the den- 
sity operator, from Eqs. (|l|) and (E^) it is apparent that 
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